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SECTION  1 
INTRODUCTION 


1.1  THE  FOIL  SWITCH  PROBLEM. 

The  Eroding  Foil  Switch  (EFS)  model  describes  foil  motion  in  a 
foil  switch.  The  foil  switch  is  a  promising  concept  for  sharpening  the 
risetime  and  thus  increasing  the  power  of  capacitor  banks  used  in  pulsed 
power  applications.1  This  report  describes  both  the  foil  switch  device  and 
experimental  results  obtained  with  it.  The  EFS  model  will  be  derived  and 
shown  to  provide  a  reasonable  explanation  of  experimental  results.  The  EFS 
model  will  also  be  shown  to  suggest  methods  of  improving  switch  performance. 

1.2  THE  FOIL  SWITCH. 

A  foil  switch  consists  of  a  coaxial  transmission  line  with  a  thin 
annular  foil  mounted  between  the  inner  and  outer  electrodes  (Figure  1). 
When  a  voltage  is  applied  by  a  capacitor  bank  to  the  electrodes,  current 
flows  in  the  foil,  and  a  Lorentz  force  pushes  the  foil  down  the  transmission 
line.  A  gap  in  the  outer  electrode  of  the  line  allows  rapid  switching  of 
current  when  the  foil  passes.  The  load  that  the  current  is  to  be  switched 
into  is  connected  across  the  gap,  so  that  the  passage  of  the  foil  across  the 
gap  will  suddenly  introduce  the  load  into  the  circuit. 

Three  important  effects  must  be  considered  when  designing  a  foil 

switch: 


1)  The  foil  should  reach  the  gap  at  high  velocity  so  that  the 
switching  time, 

it  =  d/Vf  ,  (1.1) 

where  d  is  the  gap  width  and  \f  is  the  foil  velocity  at  the 
gap,  is  as  short  as  possible. 
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2)  The  gap  must  be  kept  clear  of  both  foil  fragments  and  plasma, 
so  that  It  will  not  "short  out,"  allowing  current  to  bypass 
the  load.  This  requirement  of  a  "clean  gap"  means  that  the 
foil  must  be  shaped  so  that  it  accelerates  in  one  piece  rather 
than  fragmenting  due  to  unequal  magnetic  forces.  It  also 
requires  that  the  foil  not  be  too  light,  or  it  will  vaporize 
due  to  resistive  heating  before  it  reaches  the  gap,  filling 
the  gap  with  vapor. 

3)  The  foil  motion  must  be  calculated  to  allow  the  foil  to  reach 

the  gap  at  the  time  of  maximum  current,  that  is,  the  quarter 

cycle  time  of  the  bank.  The  inductance  of  the  switch  is 

usually  chosen  to  be  much  less  than  the  inductance  of  the  load 
or  the  driving  circuit,  so  that  the  bank  waveform  is 
unaffected  by  the  foil  motion. 

1.3  REPORT  OVERVIEW  AND  MOTIVATION. 

The  optimum  design  of  a  foil  switch  requires  understanding  of  the 
foil  motion  and  heating.  The  same  currents  that  accelerate  the  foil  will 

also  heat  it.  Therefore,  the  requirements  of  high  foil  velocity  and  no  foil 

vaporization  must  be  balanced  against  one  another.  The  foil  equation  of 
motion  also  determines  the  choice  of  gap  placement  in  the  transmission 
line,  so  that  the  foil  can  arrive  at  the  gap  at  the  time  of  maximum  bank 
current.  Present  models1,2  of  foil  motion  assume  uniform  heating  and 
constant  mass  for  the  foil.  These  models,  however,  appear  inadequate  when 
their  predictions  are  compared  with  experimental  data.  The  foils  appear  to 
reach  a  much  higher  temperature  than  is  predicted  using  a  uniform  current 
model  and  to  arrive  at  the  gap  earlier  than  expected  using  constant  mass. 

\  new  model  for  foil  motion,  which  is  far  more  consistent  with 
experimental  data,  has  been  formulated.  This  model  is  called  the  Eroding 
Foil  Switch  (EFS).  The  EFS  model  assumes  two  phases  of  foil  motion:  an 
initial  "rocket"  phase  and  a  secondary  "snowplow"  phase.  The  EFS  model 
assumes  that,  initially,  currents  in  an  accelerated  foil  flow  in  a  thin 
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layer  rather  than  uniformly.  This  effect,  due  to  nonlinear  magnetic 
diffusion,  results  in  ablation  or  erosion  in  thin  layers  of  the  foil  as  it 
accelerates.  This  is  called  the  "rocket"  phase  because  the  mass  of  the  foil 
becomes  smaller  as  it  accelerates.  The  mass  decreases  at  a  rate 
proportional  to  the  resistive  heating.  Therefore,  the  foils  in  this  model 
arrive  at  the  switching  gap  much  sooner  than  is  predicted  by  constant  mass 
models.  The  erosion  of  the  foil  can  lead  to  its  complete  disassembly.  The 
magnetic  field  pushing  the  foil  then  diffuses  through  the  exploded  foil  as  a 
shock  wave.  This  is  called  the  "snowplow"  phase. 

In  the  body  of  this  report,  we  will  first  discuss  the  experimental 
data  that  first  suggested  the  EFS  model.  A  derivation  of  the  equations  of 
motion  of  two  phases  of  the  EFS  model  will  then  be  shown.  Finally,  the  EFS 
model  will  be  applied  to  the  results  of  a  foil  switch  experiment. 


SECTION  2 

SWITCH  EXPERIfCHTS 

2.1  THE  FOIL  SWITCH  EXPERI^NT. 


!WWT' 


The  EPS  model  is  based  on  the  results  of  two  experiments.  The  most 
important  of  these  experiments  is  the  foil  switch.  The  foil  switch 
experiment1  consisted  of  a  coaxial  transmission  line  of  two  electrodes  (Figure 
1).  The  inner  electrode  had  an  0.8  cm  radius  and  an  outer  electrode  which 
formed  a  sleeve  with  an  inner  surface  radius  of  1.1  cm.  The  annular  cavity 
thus  formed  had  a  height  of  0.3  cm.  In  this  cavity,  an  annular  foil  of  0.3  cm 
height  and  0.07  cm  average  thickness  was  mounted.  A  gap  0.3  cm  wide  was  cut 
in  the  outer  electrode  beginning  1.0  cm  from  the  mounting  point  of  the  foil. 

Because  a  variation  of  magnetic  force  occurred  on  the  foil  due  to 

the  radial  variation  of  the  magnetic  field,  the  foil  was  "mass  graded"  in 

radius  so  that  the  mass  per  unit  area,  y  =  poo,  where  u  is  the  mass  per 

unit  area,  p  is  the  foil  material  mass  density,  and  6  is  the  foil  thick- 
o 

ness,  was  proportional  to  the  magnetic  force  per  unit  area  on  the  foil, 
B2/8ir  (in  cgs  units),  where  B  is  the  magnetic  field  strength.  Assuming  the 
mass  per  unit  area,  u,  was  constant  in  time,  the  acceleration  of  each  part  of 
the  foil  was  given, 

32tz  =  B2(t,r)/(8iry(r) )  .  (2.1) 


Where  z  is  the  distance  in  the  direction  of  foil  motion,  the  z  axis,  and  B 
varies  as  1/r,  due  to  the  coaxial  geometry: 


where  I  is  the  current,  c  is  the  speed  of  light,  and  r  is  the  radius.  The 
foil  mass  per  unit  area  (or  thickness,  6,)  had  to  vary  by  approximately  1/r2 
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32  z  =  A(t)  ,  (2.3) 

that  Is,  an  acceleration  that  is  constant  in  radius  and  moves  each  part  of 
the  foil  equally.  This  mass  grading  is  necessary  because  of  the  high 
accelerations  experienced  by  the  foil.  These  accelerations  would  create 
stresses  exceeding  the  foil  yield  strength  if  they  were  not  distributed 
equally.  The  foil  mass  and  distance  to  the  gap  were  also  chosen  carefully, 
assuming  that  the  foil  mass  would  remain  constant. 

The  foil  was  heated  by  the  electric  current  flowing  in  it  and 
would  melt  and  even  vaporize  if  accelerated  long  enough.  The  experimental 
parameters  were  chosen  so  that  the  foil  would  only  become  molten  (660°C) 
when  it  reached  the  gap.  These  calculations  assumed,  however,  that  the  foil 
would  heat  uniformly. 

The  actual  calculation  of  foil  parameters,  accelerations,  and 
heating  rates  will  only  be  summarized  in  this  report  but  are  calculated  in 
detail  in  a  companion  report,  "Physical  Properties  and  Numerical  Models  of 
Foils  in  High  Current  Switches."3 

The  foil  equation  of  motion,  assuming  the  foil  retained  its 
original  mass,  is  found  (using  the  characteristic  time,  t  =  10"5sec>  for 
convenience) 


3tV  =  3.7cm/ x2( sin 2wt) 


(2.4a) 


where  oj  =  2.8 t-1 


(2.4b) 


The  functions  for  z  and  V  can  be  integrated  and  become 
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V  =  1.85cm/x2[t  -  —  sin(2uit)} 
2oj 


(2.5a) 
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z  =  1.85cm/i2(t2/2  +  -  sin2(u:t))  . 

2*2 


(2.5b) 
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Similarly,  the  function  for  foil  temperature  is  found  to  be 

T(t)  =  1.5  x  107  C°/sec  (t  -  1-  sin(2ujt))  .  (2.6) 

2bl 

The  results  of  the  experiment  were  as  follows.  The  foil 
apparently  reached  the  gap  in  t  =  0.7 t  =  7x10“^sec  and  achieved  a  velocity 
of  V  =  3.0cm/T  =  3x105cm/sec.  The  foil  reached  the  gap  sooner  than  expected 
indicating  50%  of  its  original  mass  may  have  been  lost.  The  switch  signal 
was  a  linear  rise  over  0.1  t  (10“^  sec)  (Figure  2),  rather  than  a  sharp 
rise  consistent  with  a  concentrated  source  of  current.  Temperature 
diagnostics  showed  that  the  foil  had  reached  a  4000°C  temperature  rather 
than  the  expected  600°C.  Therefore,  the  foil  appeared  to  have  lost  mass  as 
it  accelerated  and  also  to  have  heated  more  strongly  than  expected. 

2.2  THE  MAGNETIC  GATE  EXPERIMENT. 

The  magnetic  gate  experiment2  was  a  switching  experiment  of 
different  geometry  from  the  foil  switch.  It  yielded  much  data  on  foil 
motion,  because  the  foil  trajectory  was  more  thoroughly  diagnosed  than  in 
the  foil  switch  experiment. 

The  magnetic  gate  experiment  differed  from  the  foil  switch  in  many 
respects,  chief  of  which  was  having  several  gaps  instead  of  one.  This  pro¬ 
vided  data  on  the  location  of  the  foil  at  several  points  in  time.  The 
magnetic  gate  experiment  also  differed  from  the  foil  switch  in  that  it  was  a 
linear  transmission  line  rather  than  a  coaxial  system  (Figure  3),  and  the 
foil  being  accelerated  consisted  of  a  thinned  region  of  a  larger  metal  gate, 
rather  than  a  simple  thin  foil.  The  gate  thickness  was  0.2  cm,  much  thicker 
than  the  foil  in  the  foil  switch.  The  foil  was  driven  by  a  capacitor  bank 
delivering  2.5  MA  peak  current  with  a  quarter  cycle  time  of  4.8  micro¬ 
seconds. 


The  magnetic  gate  was  designed  to  open  by  stretching  out  the 
thinned  region  when  the  magnetic  pressure  exceeded  the  yield  strength  of  the 
gate.  The  thick  portions  of  the  gate  would  then  bend  down  as  hinges  and 
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Figure  2.  Plot  of  normalized  switched  current 
(switched  current  divided  by  bank 
current)  for  the  Los  Alamos  foil 
switch  experiment. 


8 


'“.'••.v.v.-.-.  .*  -v 


allow  the  thin  region  to  accelerate  as  a  foil.  The  resulting  current 
carrying  foil  switched  the  current  through  the  gaps  with  very  high 
efficiency,  switching  almost  100%  of  the  bank  current  (Figure  4).  However, 
the  foil  moved  as  If  it  had  lost  95%  of  its  original  mass  (Figure  5);  the 
foil  moved  down  the  transmission  line  much  more  rapidly  than  anticipated.  In 
particular,  the  last  data  point  on  the  plot  of  foil  position  versus  time 
indicates  it  is  losing  mass  as  it  accelerates.  This  is  because  the  last  data 
point  is  above  the  curve  fitting  the  first  two  points  and  thus  indicating 
that  it  is  moving  faster  than  expected. 

2.3  THE  INTERPRETATION  OF  EXPERIMENTAL  DATA. 

Despite  the  differences  between  the  foil  switch  and  magnetic  gate 
experiments,  similar  phenomena  occurred  in  each.  In  both  experiments,  large 
mass  losses  were  observed  in  the  accelerated  foils.  In  addition,  foil 
heating  was  observed  greatly  in  excess  of  what  was  expected  from  models  of 
uniform  current  flow  in  the  foils. 

In  the  foil  switch  experiment,  the  kinematic  results  were 
consistent  with  roughly  a  50%  initial  mass  loss  from  the  foil.  However,  an 
examination  of  the  apparatus  after  the  experiments  revealed  a  continuous  film 
of  foil  residue  down  the  line,  rather  than  any  concentration  of  debris  at  the 
foil  starting  point  that  would  indicate  an  abrupt  initial  mass  loss.4  The 
experimenters  suggested  viscous  entrainment  of  foil  material  as  a  mechanism 
for  mass  loss. 

The  switch  signal  itself  was  not  what  one  would  expect  from  a 
concentrated  current  flowing  in  a  foil.  However,  the  switch  current  signal 
could  be  explained  by  a  simple  physical  model.  The  model  consisted  of  a 
region  of  uniform  current  density  1  cm  wide  that  moved  at  uniform  speed 
across  the  gap.  The  current  density  acted  like  a  commutator  brush  to 
transfer  current  from  one  side  of  the  gap  to  another.  The  switch  current 
through  the  gap  was  then 

Is  =  30x  ,  (2.7) 
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Figure  5.  A  plot  of  foil  position  versus  time  from 
the  Magnetic  Gate  Experiment.  (1)  ex¬ 
periment  data  (squares),  (2)  F=ma  calcu¬ 
lations  using  the  entire  gate  mass 

(dotted  curve),  (3)  F=ma  calculations 
using  5%  of  the  original  gate  mass  (solid 


where  30  is  the  current  density  in  megaamps  per  cm;  x  is  defined  as 

x  =  V0t  -  d  ,  (2.8) 

where  V0  is  the  uniform  velocity,  and  d  is  the  location  of  the  far  side  of 
the  gap.  This  expression  obviously  gives  a  linear  rising  current  as  is  seen 
in  the  experiment.  However,  the  uniform  current  is  spread  over  a  large 
(-1  cm)  area,  which  was  much  larger  than  the  foil  mean  thickness  at  its 
starting  point.  If  the  model  of  a  diffuse  uniform  current  was  correct,  then 
the  foil  was  no  longer  a  dense  object,  but  a  vapor.  This  was  consistent 
with  optical  pyrometry  signals  from  the  foil  which  indicated  a  4000°C 
temperature.  The  foil  had  evidently  vaporized.  However,  using  the  model 
that  the  foil  would  carry  a  uniform  current  and  thus  heat  uniformly,  the 
foil  mass  had  been  chosen  to  ensure  that  it  would  not  reach  such  high 

temperatures. 

In  the  magnetic  gate  experiment,  a  large  mass  loss  was  observed  as 
in  the  foil  switch  experiment.  In  addition,  the  foil  was  observed  to  become 
incandescent,  indicating  high  temperature,  though  precise  temperature 
measurements  were  not  made. 

Therefore,  the  two  experiments  indicate  that  models  of  foil 
motion,  which  assume  uniform  current  density  and  constant  foil  mass,  are 

inadequate.  A  new  model  for  foil  motion  allowing  mass  loss  and  rapid 

heating  is  needed.  Such  a  model  will  now  be  derived  based  on  a  careful 

study  of  the  flow  of  current  in  strongly  accelerated  foils.  The  effects  of 
nonlinear  magnetic  diffusion  and  strong  acceleration  will  be  found  to  result 
in  both  mass  loss  and  rapid  heating  of  foils. 
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SECTION  3 

THE  PHYSICS  OF  THE  EFS  MODEL 


In  the  previous  section,  we  discussed  the  results  of  experiments 
on  strongly  accelerated  foils.  In  both  of  these  experiments,  the  mass  of 
the  foils  appeared  to  decrease  during  acceleration,  and  the  foils  became 
much  hotter  than  was  expected.  Both  of  these  effects  can  be  explained  if 
current  did  not  flow  uniformly  in  the  foil  but  instead  concentrated  in  a 
small  region  which  then  moved  through  the  foil. 

If  the  current  concentrated  in  a  layer,  thin  compared  to  the 
thickness  of  the  foil,  the  current  would  strongly  heat  the  foil  locally. 
Also,  the  force  pushing  the  foil  would  be  concentrated  in  this  layer,  and  as 
the  layer  moved  forward,  it  would  "leave  behind"  heated  foil  material.  The 
foil  material  behind  the  concentrated  current  would  no  longer  be  accelerated 
but  would  instead  be  under  tension,  and  if  accelerations  were  strong,  would 
be  tom  off  the  foil  in  the  form  of  small  droplets.  The  foil  mass  would 
then  decrease  as  the  foil  was  accelerated. 

In  this  section,  the  phenomenon  of  nonlinear  magnetic  diffusion 
will  be  shown  to  lead  to  concentration  of  current  in  a  thin  foil.  The 
effects  of  this  current  concentration  will  be  shown  to  lead  to  mass  erosion 
in  strongly  accelerated  foils.  The  effect  of  mass  erosion  on  the  foil 
equation  of  motion  will  then  be  discussed.  This  erosion  can  consume  the 
foil  completely.  The  foil  material  will  then  explode  into  vapor  with  the 
magnetic  field  moving  through  it  as  a  shock  wave.  The  EFS  model  is  then 
seen  to  arise  easily  from  these  concepts.  The  EFS  model  can  be  said  to  be 
the  kinematics  of  light  foils  undergoing  nonlinear  magnetic  diffusion. 

3.1  NONLINEAR  MAGNETIC  DIFFUSION:  NUMERICAL  AM)  APPROXIMATE  TREATMENTS. 

When  a  magnetic  field  Is  applied  to  any  conductor,  all  current 
initially  is  in  a  thin  layer  on  the  conductor  surface.  As  time  progresses, 
the  current  layer  broadens  and  fills  the  conductor  interior.  The 
concentration  of  current  in  a  thin  layer  will  result  in  a  strong  ohmic 
heating  rate, 

3tQ  =  nO2  ,  (3.1) 
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where  Q  is  the  thermal  energy  per  unit  volume  in  ergs/cm  ,  3  is  the  electric 
current  density, 

3  --  1  3  B  ,  (3.2) 

4  it  z 

from  Amperes  law  in  our  geometry,  and  n  is  the  resistivity  in  seconds  (cgs 
units).  In  general,  n  is  a  function  of  Q, 

n  —  +  CjQ  +  C2Q  ...  (3.3) 

The  constants  and  C2  were  chosen  to  fit  a  curve  in  reference  5.  They 
were  given  the  values 

CA  =  4.8x10~2i  sec/bar  and  C2  =  3.9x10" 24  sec/bar2 
by  the  authors  for  the  purposes  of  numerical  calculations. 

Ignoring  radial  variations,  we  can  write  using  Faraday's,  Ohm's, 
and  Ampere’s  laws, 

c2 

3tB  =  —  3z(n3zB)  .  (3.4) 

4u 

Realistic  resistivity  models5  were  employed  for  Equations  3.1, 
3.3,  and  3.4,  which  were  solved  as  a  coupled  system.  Results  for  a  magnetic 
field-foil  system  similar  to  the  experiment  (Figure  6)  showed  the  current 
forming  a  layer  and  propagating  as  a  wave  through  the  system.  Since  3xB 
forces  would  be  concentrated  in  this  layer,  the  metal  left  behind  by  the 
current  layer  passage  would  no  longer  be  pushed  by  the  magnetic  pressure  and 
thus  be  left  behind  (assuming  that  being  heated,  it  would  lose  its  material 
strength  under  these  conditions  and  slough  or  spall  off  in  droplets). 

Analytically,  the  advance  of  the  diffusion  wave  can  be  understood 
by  studying  the  set  of  Equations  3.1  -  3.4  and  making  some  physically 
reasonable  approximations.  A  more  general  approximate  treatment  of  this 
problem  is  given  in  the  next  section  of  this  report;  also,  an  extensive 
discussion  of  this  effect  is  given  in  Reference  6. 


The  current  begins  as  a  thin  layer  decaying  exponentially  into  the 
unheated  metal  with  a  characteristic  length  of  a  "skin  depth," 


3  =  30ei 


2c  A/2 

•  « *  (4^) ' 

4ttoj 


(3.5) 


where  u  is  the  field  frequency  or  inverse  rise  time.  As  the  metal  heats  due 
to  the  skin  current,  it  will  become  more  resistive,  and  the  current  will 
also  diffuse  more  deeply  into  the  metal.  The  current  will  now  decay  behind 
the  layer  of  peak  current  density,  because  the  resistivity  will  be  high. 
The  current  will  thus  concentrate  where  resistivity  is  just  beginning  to 
rise  in  the  metal.  In  this  region,  the  resistivity  will  be  nearly  linear  in 
internal  energy,4  so  that  Equation  3.3  reduces  to: 


n  *  n  +  C|Q 
o 


(3.6) 


We  can  thus  write  for  Equation  3.1: 


n  =  C|T"i3  • 


(3.7) 


We  also  consider  this  effect  of  a  peak  in  current  density  in  the  region  of 
small  increase  in  resistivity:  in  the  region  of  interest,  the  magnetic 

field  will  increase  strongly,  but  then  assume  nearly  the  value  it  has  on  the 
boundary.  In  the  region  of  interest,  where  the  diffusion  wave  is  actually 
advancing,  we  have  a  peak  of  current, 


341=0  , 


(3.8a) 


which  means: 


3£B  =  0 


(3.8b) 


We  can  then  write  for  Equation  3.4 


3+-B  -  —  3^ ri  3ZB  • 
4ir 


(3.9) 


This  equation  describes  a  convection  of  B  across  the  system  in  a 
self- similar  manner  at  a  speed: 

jl 


VB  =  ~  3Zn 


(3.10) 


We  now  can  solve  for  this  speed  by  solving  Equations  3.7  and  3.9  together 

c2C,  ? 

3tn  =  3Bn3tB  =  - -  ( 3ZB) 2  ,  (3.11a) 

(4tt)2 


assuming  n  and  Q  are  functions  of  B: 

3Bti3zti3zB  =  cir|(3TB)2  > 

which  simplifies  to: 

(3zn)2  «  £iH  (3ZB)2 
4ir 

Integration  in  z  over  the  region  of  linear  behavior  gives 

n  q2 


(3.11b) 


(3.11c) 


(3.12) 


This  relation  can  be  rewritten: 


AQ  =  AB2/8ir  . 


energy  densities  couple  at  just  the  right  rate  to  balance  each  other.  That 
is,  energy  density  changes  are  shared  equally.  This  follows  as  the  result 
of  the  approximations  in  the  region  of  dynamics  given  by  Equations  3.6  and 
3.7.  The  Equipartition  of  energies  is  an  interesting  physical  result  that 
may  be  more  general  than  the  approximations  used  to  find  it.  When 


k  W  J  »  .  ■«  *  *  m  ’  .  *.•.*,  *  .  *  »  ■  k  -  -  -  »  -  -  -  -  •  -  - 

>  -  .  •  .n  .  •  ,  *  ,  * V  .*  •*>  *  k  -  .  • .  - 


>  .v  .% 


(3.13) 

O 

that  in 

* » • , 

thermal 

Equipartition  is  used  itself  as  an  initial  assumption,  it  is  found  to  allow 
great  simplification  of  the  nonlinear  diffusion  problem.  This  approach  will 
be  discussed  in  the  next  section. 

3.2  NONLINEAR  MAGNETIC  DIFFUSION:  THE  EQUIPARTITION  APPROXIMATION. 

We  attempt  to  find  a  general  solution  for  the  nonlinear  magnetic 
diffusion  problem  that  is  insensitive  to  the  exact  dependence  of 
resistivity,  n,  on  the  thermal  energy  density,  Q,  providing  that  n  is 
monotonically  increasing  with  Q. 

We  find  this  general  solution  for  magnetic  diffusion  by  making  an 
approximation  that  is  consistent  with  the  case  of  a  linear  resistivity 
function.5  The  approximation  used  to  find  a  general  solution  can  be  termed 
an  "Equipartition  Approximation;"  that  is,  energy  in  the  system  should  tend 
to  divide  itself  equally  into  magnetic  and  thermal  energies.  Therefore,  we 
should  have  generally: 

B2/8n  =  Q  .  (3.14) 

This  can  be  seen  to  be  a  reasonable  assumption  by  Poyntings  Theorem  in  the 
limit  of  negligible  foil  motion: 

3.(B2/8ir  +  Q)  =  £_  3  (EB)  ,  (3.15) 
Z  4x  Z 

which  indicates  that  both  magnetic  and  thermal  energy  densities  can  be 
lumped  together  with  one  source  term.  The  suggestion  that  the  source  term 
would  supply  both  energy  densities  equally  seems  an  obvious  one. 

Keeping  this  approximation  in  mind,  we  can  write  both  the  magnetic 
diffusion  equation  and  the  heating  equation  in  terms  of  the  electric  field, 
E,  through  Ohm's  law, 

E  =  n3  (3.16) 
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and  Ampere's  law, 


cr> 

E  =  —  3  B  .  (3.17) 

4  it  z 


Writing  the  magnetic  equation,  Equation  3.4,  as: 
3tB  =  -c  3ZE  , 


(3.18) 


and  heating  equation,  Equation  3.1,  as: 

3tQ  =  £_  E32B  , 

4it 


(3.19) 


we  can  use  our  Equipartition  Approximation  to  eliminate  Q  and  use  Equation 
3.18  to  write  for  Equation  3.19: 

-£  C3ZE  =  -£  E3ZB  .  (3.20) 

4x  4  it 


We  can  thus  find  the  relation: 

InE  =  InB  +  lnL(t)  (3.21) 

which  becomes 

E  =  L(t)B  ,  (3.22) 

where  L(t)  is  a  function  of  time  only.  Because  E  contains  the  resistivity, 
the  exact  form  of  the  resistivity  is  not  important  as  long  as  n  rises  with 
Q. 


We  can  then  write  for  Equation  3.18: 
3tB  =  -L(t)3zB  . 


(3.23) 
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The  Equipartition  Approximation  has  therefore  allowed  us  to  take  the  general 
nonlinear  diffusion  problem  and  reduce  it  to  a  linear,  first  order 
differential  equation  for  B.  Equation  3.22  is  merely  the  equation  for  the 
translation  of  B  at  a  speed  l(t).  Defining 

z0(t)  =  Jo  L(t')df  ,  (3.24) 


we  can  write 

B  =  B0F(z0(t)  -  z) 


where  F  is  defined  through 

F(zo(t»  =  Sill  L  ,  o 

Dq  1 

consistent  with  our  less  general  treatment, 

L(t)  =  y  — (t)  . 

Sir 


We  obtain  for  a  sinusoidal  magnetic  field: 

2  (t)  - sln2’* )  . 

0  2  Sit  2ui 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


We  also  can  obtain  an  approximate  solution  following  Knoepfel5  by  assuming 
linear  resistivity  and  Equipartition  and  the  approximation  that 

B  =  B0  (t/c)1/2  (3.29a) 

with 

t  =  Vgt-z  (3.29b) 


with  Vg  a  constant. 


8*1 


It  follows  that 


Sir 


(3.30) 


where  l  is  the  characteristic  length  over  which  B  varies  in  the  region  of 
sharp  increase.  In  order  for  the  two  approximate  treatments  to  agree 
asymptotically,  we  must  have  for  y  in  Equation  3.27: 
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c2Ci 

4iU 


(3.31) 


The  diffusion  wave  speed  depends  on  the  magnetic  field  strength 
squared  and  thus  the  resistive  heating  rate.  The  diffusion  wave  velocity  will 
then  be  approximately, 


VB 


32  n 


c2CxB2(t) 


(3.32) 


where  we  assume  i  =  c/2. 

Using  c  =  0.04  cm,  CA  »  4.8  x  10"21  sec/bar,  we  obtain 


VB 


8cm  B02 

sec- bar  8tt 


(3.33) 


where  B0  is  the  magnetic  field  at  the  boundary  and  where  a  bar  is  10 6 
ergs/cm3.  This  compares  very  favorably  with  the  value  obtained  from 
numerical  solution  of  Equations  3.1  through  3.4. 


in  02 

10  cm  o 

B  sec-bar  8ir 


(3.34) 


These  values  will  vary  weakly  with  <d  through  c,  however,  given  that  id  varies 
only  slightly  between  experiments  c  =  .04cm  can  be  considered  a  constant. 

3.3  MASS  EROSION  AND  MAGNETIC  ACCELERATION:  THE  "ROCKET"  EQUATION. 


The  nonlinear  magnetic  diffusion  will  lead  to  mass  erosion  in 
strongly  accelerated  foils  if  the  current  layer  formed  is  thin  compared  to  the 
foil  thickness.  This  must  occur  because  the  "magnetic  pressure"  or  Lorentz 
force  will  push  only  that  foil  that  is  in  front  of  the  current  layer.  The 
foil  material  already  heated,  resistive,  and  passed  by  in  the  advance  of  the 
current  layer  will  be  torn  away.  This  effect  will  only  occur  in  thin  foils 
undergoing  high  acceleration.  Massive  blocks  of  metal  may  only  suffer 
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sloughing  off  of  metal  on  long-time  scales.  The  effect  may  also  disappear 
in  very  thin  foils  where  the  current  layer  dimension  is  comparable  to  the 
thickness. 


In  a  strongly  accelerated  foil  the  time  derivitive  of  u,  the  foil 
mass  per  unit  area,  is  then: 

3tu  =  -P0VB  (3.35) 

where  pQ  is  the  mass  density  of  aluminum,  2.7  gm/crn  . 

We  can  define  K: 

,  B2(t) 

3t  u  =  -K  — .  (3.36) 


K  =  P0T 

and  evaluate  it,  approximately,  using  the  value  for  Vg  obtained 

analytically: 

K=  — -ll9- -  .  (3.37) 

cm2  sec-bar 

What  then  is  the  effect  of  mass  erosion  on  foil  equations  of 
motion?  The  answer  is  found  by  solving  the  equations  for  motion  for  each 
part  of  the  foil.  For  variable  mass,  Newton's  law  is  written: 

3t(uV)  =  B2/8tr  .  (3.39) 

We  integrate  Equation  3.38  in  time: 

uV  =  W(t)  ,  (3.39) 

where  W(t)  is  an  "action  density," 

W(t)  =  /B2/8tt  dt  .  (3.40) 
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From  Equation  3.36a  we  have  for  u, 


3-£  -KB  /8*n 


(3.41) 


we  integrate  both  sides: 


U  =  -KW(t)  +  u0  i 


(3.42) 


where  Uq  is  the  initial  mass  per  unit  area.  We  can  write 


u  =  u  11  -  K  —  —  1 
*  O  V  Wo 


We  can  then  write  for  V: 

W(t) 


V  = 


,  W(t). 

“of1  -  K  — ) 


(3.43) 


(3.W) 


$ 

$ 


This  is  the  equation  for  velocity  in  a  system  where  mass  erosion 
occurs  at  a  rate  proportional  to  magnetic  pressure.  This  equation  is  quite 
general  and  allows  any  time  variation  of  the  magnetic  field  to  be  used.  The 
equation  of  motion  is  termed  the  rocket  equation,  because,  like  a  rocket, 
the  mass  decreases  as  the  foil  accelerates. 

3.4  HASS  EROSION  COUPLED  WITH  MASS  GRADING. 


■  V 
*  * 
•  V 


One  important  effect  that  must  be  considered  is  the  effect  of  mass 
erosion  coupled  to  the  variation  of  magnetic  pressure  across  the  foil 
surface.  Will  mass  erosion  destroy  the  effects  of  mass  grading? 

The  foil  switch  is  based  on  one-dimensional  motion;  that  is,  it 
assumes  the  foil  will  move  essentially  as  a  flat  sheet  rather  than 
distorting  or  bursting.  In  the  geometry  of  the  experiment,  the  magnetic 
field  varies  as  1/r,  and  thus  the  magnetic  pressure  varies  as  1/r2.  The 
foil  was  also  mass  graded  so  that 


1C. 


& 


ft 


=  B2/(8iTy)  =  A(t) 


(3.45) 


* 

s.\n 

>5a 
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therefore,  the  whole  foil  should  accelerate  as  a  unit.  This  was  confirmed 
by  simulations  of  the  experiment  performed  using  the  code  AXFOIL,  which 
simulated  the  foil  in  two-dimensions  as  a  collection  of  point  masses. 
However,  the  mass  grading  presumed  that  y  would  be  constant  in  time.  If  m 
varied  in  time  by  Equation  3.41a,  would  this  cause  even  a  mass-graded  foil 
to  be  unstable  and  distort  in  two-dimensions? 


We  use  Vc,  the  velocity  for  a  constant  y  (Vc  =  W(t)/uo),  and  write 


V(t) 


Vc(t) 

(1  -  KVc(t)) 


(3.46) 


For  the  case  of  a  mass-graded  foil, 

3rVc(t)  =  0  ,  (3.47) 

Vc  is  a  constant  in  radius;  that  is,  all  the  parts  of  the  foil  accelerate 
together,  and  accordingly, 

3rV(t)  =  0  .  (3.48) 

Therefore,  mass  erosion  at  a  rate  proporational  to  acceleration,  allows  mass 
grading  and  thus  one  dimensional  motion  to  be  preserved. 

3.5  APPROXIMATIONS  AM)  PROBLEMS  ASSOCIATED  WITH  MASS  EROSION. 

We  have  seen  in  this  section  that  nonlinear  magnetic  diffusion  and 
strong  acceleration  can  combine  to  give  a  new  picture  of  foil  motion:  a 
model  where  mass  disappears  and  acceleration  is  much  more  rapid  than 
expected.  We  have  found  the  equations  of  motion  for  mass  erosion  show  that 
this  is  consistent  with  mass  grading  to  achieve  one-dimensional  motion.  It 
must  be  pointed  out,  however,  that  this  last  statement  is  only  approximately 
true.  The  scale  lengths  for  magnetic  diffusion  will  not  change  with  mass 
grading,  so  that  the  implied  approximation  used  to  describe  mass  erosion, 
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that  foil  thickness  is  much  greater  than  the  diffusion  scale  lengths,  may 
not  be  valid  in  the  thinnest  regions  of  a  mass-graded  foil. 

We  have  another  problem  with  the  mass  erosion  effect  that  has  been 
derived.  The  foil  in  the  LANL  foil  switch  experiment  will  erode  completely 
before  it  reaches  the  gap  if  the  erosion  rates  we  have  derived  are  correct. 
Also,  the  experimental  data  indicated  a  diffuse  vapor  passing  through  the 
gap,  not  a  rapidly  moving  foil. 

In  order  to  have  a  complete  and  self-consistent  picture  to 
interpret  the  foil  experimental  results,  a  picture  of  events  after  mass 
erosion  has  destroyed  the  foil  must  be  added.  The  model  will  be  discussed 
in  detail  in  the  next  section. 

3.6  FOIL  MOTION  AFTER  EROSION:  THE  SNOWPLOW  EQUATION. 

As  we  have  seen,  foil  erosion  can  lead  to  a  rocket- like 
acceleration.  However,  the  mass  erosion  rate  we  have  derived  indicates  that 
foils  in  many  experiments  will  erode  completely.  According  to  Equation 
3.44,  this  will  give  an  infinite  velocity  at  y=0.  This,  obviously,  does  not 
occur.  What  will  occur  will  be  a  "saturation"  of  the  velocity  at  some  high 
value  with  acceleration  ending. 

The  exact  process  of  this  saturation  Is  difficult  to  characterize 
in  detail,  however,  it  can  be  described  qualitatively. 

Before  the  foil  mass  becomes  zero  (y=0)  the  last  remaining  layer 
of  foil  will  heat  rapidly  and  explode.  This  will  occur  because  the 
nonlinear  diffusion  process,  while  it  can  occur,  actually  prevents  any 
portion  of  the  foil  from  having  a  thermal  energy  density  above  B2/8it,  i.e., 
Q  *  B  / 8 it.  Once  the  advancing  nonlinear  wave  reaches  the  front  surface  of 
the  foil,  however,  the  current  will  have  no  further  place  to  diffuse.  The 
remaining  layer  of  foil  must  then  heat  until  it  explodes  into  vapor.  This 
is  what  seems  to  occur  in  the  LANL  Foil  Switch  Experiment. 
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We  will  assume  the  vaporized  foil  will  explode  forward  ahead  of 
the  magnetic  field  and  fill  the  gap  uniformly. 

The  vapor  will  be  partly  ionized  and  thus  be  a  conductor.  The 
magnetic  field  will  advance  through  it  like  a  shock  wave.  We  can  describe 
this  mathematically  by  using  Newton's  law: 

V3tu  +  y3tV  =  B2/8it  .  (3.48) 

3t(wV)  =  B2/8w  (3.49) 

In  this  case  we  assume  that  material  will  be  swept  up  by  the  shock 
and  then  diffuse  through  it  after  being  brought  up  to  speed.  Therefore,  we 
have  the  approximation 

3ty  »  y3tV  .  (3.50) 

We  then  have 


V3ty  =  B2/8it  (3.51) 

with 

3fcU  =  PvV  ,  (3.52) 

where  pv  =  the  mass  density  of  the  exploded  foil  cloud.  We  can  then  solve 
for  V  using  Equations  3.51  and  3.52: 

-  R2  T1/2 

vs  =  !— L_  ]  .  (3.53) 


where  Vs  is  "snowplow"  velocity.  The  velocity  "saturates".  No  longer  is 
acceleration  proportional  to  B  .  Instead,  velocity  has  become  proportional 
to  magnetic  field  strength. 


27 


3.7  DISCUSSION 


In  this  section  we  have  derived  the  equations  of  motion  for  light  foils 
undergoing  nonlinear  magnetic  diffusion  and  high  acceleration.  We  saw  how 
this  could  lead  to  mass  erosion  and  rapid  acceleration  of  the  foil.  This 
motion  was  described  by  the  rocket  equation.  It  was  also  found  that  foils 
could  erode  completely  and  explode  during  acceleration.  The  magnetic  field 
then  moving  as  a  shock  wave  through  the  exploded  material.  This  was 
described  by  the  snowplow  equation.  In  the  next  section  the  equations  of 
motion  will  be  used  to  define  the  EFS  model  of  foil  motion  and  apply  the 
model  to  the  LANL  switch  experiment. 


SECTION  4 

THE  EFS  MODEL  OF  FOIL  MOTION 


In  the  previous  section,  we  discussed  the  chief  physical 
phenomenon  leading  to  the  EFS  model:  the  concentration  of  current  in  the 
foil  into  a  thin  sheet  rather  than  its  uniform  spread  over  the  foil 
thickness.  In  this  section,  we  will  discuss  the  EFS  model  and  then  apply  it 
to  interpret  the  results  of  the  LANL  experiment. 

4.1  THE  "ROCKET"  AND  "SNOWPLOW"  PHASES  OF  EFS. 


The  EFS  model  predicts  two  phases  of  the  foil  motion  in  a  foil 
experiment.  The  first  phase  can  be  termed  the  mass  erosion  or  "rocket" 
phase,  because  the  equation  for  foil  motion  resembles  that  of  a  rocket;  for 
a  prescribed  thrust,  the  foil  accelerates  very  rapidly,  because  its  mass 
decreases  in  time  at  a  rate  proportional  to  the  thrust.  The  second  phase  is 
termed  the  "snowplow"8  phase,  because  the  foil,  being  largely  eroded, 
explodes  forward  and  becomes  a  mixture  of  plasma,  vapor,  and  liquid  droplets 
in  the  remainder  of  the  transmission  line.  The  magnetic  field  moves  through 
this  exploded  foil  material  like  a  shock  wave. 


In  the  first  phase,  the  foil  equation  of  motion  will  be  of  the 
form  discussed  in  the  previous  section: 


VR  = 


W(t) 


(4.1) 


where  K  =  21  gm/(cm  sec-bar),  and 


.  B2 

W(t)  =  /  —  (t)dt 


(4.2) 


In  the  second  phase,  a  shock  wave  or  snowplow  phase,  the  velocity 
of  the  foil  will,  in  effect,  saturate.  The  foil  is  assumed  to  explode  or 
disassemble  into  droplets  and  vapor.  To  describe  this  regime,  we  have  the 
"snowplow"  equation 


where  pv  is  the  vapor  mass  density  in  the  transmission  line.  We  thus  have 
expressions  for  the  foil  velocity  before  and  after  t£,  the  foil  explosion 
time.  To  complete  the  EFS  model,  we  assume  continuity  of  the  functions  for 
velocity  and  distance  at  t  =  tjr 

Vr  (tg;)  =  Vs(t£)  (4.4) 

These  four  equations  form  the  EFS  model  of  foil  motion.  We  can  now  apply 
this  model  to  the  LANL  foil  switch  experiment. 

K.2  INTERPRETATION  OF  THE  LANL  FOIL  SWITCH  RESULTS. 

The  LANL  switch  experiment  results  were  that  the  foil  traveled  a 
distance,  Xf  =  1.1cm,  in  a  time,  tf=7ps,  and  that  its  final  velocity  was 
approximately,  V f  =  3x10s  cm/sec.  We  can  infer  also  that  the  foil  changed 
state  during  its  acceleration,  accelerating  first  as  a  solid  object  and  then 
becoming  a  diffuse  shock  wave. 

We  wish  to  explain  these  results  in  terms  of  the  EFS  model.  Since  the 
LANL  results  are  in  terms  of  final  distance  and  velocity  we  can  conveniently 
compare  the  EFS  model  with  experiment  by  integrating  the  EFS  Equations  4.1 
and  4.3  in  time  to  obtain  the  total  distance  traveled. 


We  have  then 


X 


f 


i 

Jt^ 


L 


V  (t)  dt 


*£  W(Mdt 
-o 


u  (1-K 
o 


(4.5) 


where  Xf  is  the  total  distance  to  the  gap  and  for  the  LANL  experiment 
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(4.6) 


W(t)  =  5x103  bars  ( t- 1  / ( 2 oj)  sin  ( 2 ojt ) ) 

The  shock  velocity  is  given  by 

V  <t>  =  V  ,  (4.7) 

S  f  B(tf ) 

assuming  a  uniform  foil  debris  density.  At  this  point,  we  can  insert  our 
theoretical  values  of  the  erosion  constant,  K,  and  obtain  tE. 
Alternately,  we  can  solve  for  tE  and  K  simultaneously  and  then  compare  the 
value  of  K  obtained  experimentally  with  that  obtained  from  diffusion 
theory.  Here  we  will  solve  for  both  K  and  tE  experimentally  and  compare  K 
with  the  theoretical  value. 

We  solve  for  the  foil  break-up  time  tE  and  the  mass  erosion 
constant  K  by  solving  Equation  4.5  consistent  with  a  matching  of  shock 
velocity  to  foil  speed  at  breakup: 

V  (t.)  =  - 1 -  .  (4.8) 

w(tE) 

We  can  eliminate  K  and  solve  for  tE  using  Equation  4.9 

X  =  /  f  V  (t)  dt  ♦  f-  - — -  .  (4.9) 

f  *E  s  0  f  «o  Wo  In 

w(t)  '  W(tE)  +  vs(tE) 

When  Equation  4.19  is  solved  by  iterations,  it  gives  the  value  tE  = 
3.9x10-tsec.  The  value  of  K  obtained  is  then 

K  =  19gm/[cm^  sec-bar]  ,  (4.10) 

which  is  very  close  to  the  value  calculated  in  the  last  section, 

K  =  21gm/[cm2  sec-bar]  .  (4.11) 

The  graph  of  velocity  versus  time  for  the  foil  is  given  in  Figure  7. 
Therefore,  the  EFS  model  would  seem  to  be  a  successful  model  for  explaining 
the  LANL  foil  switch  results. 
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Figure  7 
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A  plot  of  velocity  versus  time  for  the  EFS  model  cal¬ 
culation  of  the  Los  Alamos  foil  switch  experiment. 
Note  that  the  sudden  change  in  slope  of  the  curve  at 
3.9  microseconds,  the  calculated  foil  explosion  time. 
Time  is  given  on  the  lower  axis  in  10" 5  seconds. 


SECTION  5 

DISCUSSION  AND  SUMMARY 


We  have  given  a  research  report  on  the  Eroding  Foil  Switch  model. 
We  have  demonstrated  that  the  EFS  model  is  consistent  with  the  picture  of 
nonlinear  magnetic  diffusion.  Highly  stressed  foils  may  have  very 
nonuniform  current  distributions,  rather  than  having  a  uniform  current 
density  and  heating. 

The  fact  that  current  is  localized  and  moves  through  a  foil  as  a 
thin  layer,  with  an  attendant  localization  of  3XB  forces  and  heating  in  a 
thin  layer,  means  that  the  effective  mass  of  the  foil  being  pushed  will 
decrease.  It  should  be  noted  that  this  effect  may  be  absent  from  more 
massive  pieces  of  metal  undergoing  acceleration.  If  accelerations  are 
small,  metal  behind  the  current  layer  will  remain  attached  to  the  main  mass 
of  metal,  and  mass  will  not  decrease.  However,  for  the  case  of  thin  foils, 
mass  erosion  should  be  an  an  important  effect. 

The  rate  of  mass  loss  per  unit  area  was  found  to  be 

V  =  “P0VB  (5-1} 

VQ  «  B2(t)  ,  (5.2) 

where  pQ  is  the  mass  per  unit  area.  This  is  an  entirely  reasonable 
scaling  for  mass  loss,  since  it  is  proportional  to  the  resistive  heating 
rate.  This  leads  to  "rocket"  equation  of  motion,  because  the  foil  mass 
decreases  as  it  accelerates.  The  rate  of  erosion  was  also  found  to  depend 
weakly  on  the  frequency  or  rise  time  of  the  magnetic  field,  being  faster  for 
faster  rise  times. 

A  mass  loss  rate  with  this  scaling  is  found  to  be  consistent  with 
stable  acceleration  of  mass  graded  foils.  The  EFS  model  also  suggests  that 
foils  which  would  be  expected  to  remain  solid  during  acceleration,  using 
assumptions  of  uniform  current,  may  be  vaporized  instead.  In  the  case  of 
total  foil  erosion  the  motion  of  the  foil  resembles  a  shock  wave. 


The  shock  wave  speed  is 


Vs  =  B/(81tpv)1/2  .  (5.3) 

This  is  very  similar  to  a  conventional  shock  wave  speed  with  pressure  being 
replaced  by  B  /8it,  the  magnetic  pressure.  This  is  the  "snowplow"  phase  of 
the  motion,  named  from  its  simlarity  to  an  WD  shock  wave.  This  state  may 
be  thought  of  as  a  "saturation"  of  foil  velocity,  since  acceleration  ceases 
with  the  formation  of  the  snowplow.  The  state  of  the  foil  material  in  the 
snowplow  phase  is  more  poorly  defined  than  in  the  rocket  phase.  The  only 
properties  that  can  be  inferred  are  low  mass  density  and  some  electrical 
conductivity.  More  experimental  data  defining  the  foil  state  in  this  phase 
of  motion  are  clearly  required. 

The  EFS  model  suggests  that  the  foil  switch  has  great  potential  as  a 
pulsed  power  tool.  The  model  predicts  that  a  foil  can  reach  very  high 
velocity  yet  have  low  temperature.  The  eroding  material  will  be  in  the  form 
of  liquid  or  solid  particles.  Optimum  switching  may  in  fact  occur  before 
the  foil  erodes  completely.  This  would  allow  a  high  speed  foil  to  cross  the 
gap  leaving  only  a  cloud  of  aluminum  droplets  in  the  gap. 

In  future  reports,  the  EFS  model  will  be  applied  to  help  interpret  and 
optimize  the  results  of  upcoming  experiments  which  will  be  conducted  by 
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Western  Research  Corporation.  The  EFS  model  itself  will  be  refined  as 
experimental  data  becomes  more  abundant. 
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